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During the second half of the 19th century it became a widely accepted scientific 
view that matter consists of atoms/molecules and that their motion is governed by the 
same dynamical laws as the trajectories of heavenly bodies. Of course, there was no 
knowledge on the precise form of the interaction potential. One relied on a model of 
very tiny balls which interact through elastic collisions (hard spheres), perhaps improving 
phenomenologically to a smooth interaction potential of short range. The necessity of the 
more realistic Lennard-Jones type potential was recognized only considerably later on the 
basis of, at the time, very precise thermodynamic measurements. 

In 1872 Ludwig Boltzmann started from the then accepted atomistic model and wrote 
down a nonlinear transport equation for the one-particle distribution function. In modern 
terms, he argued that at a collision the two incoming velocities are statistically indepen¬ 
dent, which leads to a closed equation for the Boltzmann / function. Such assumption 
could hold approximately only for a raribed gas. In the 6th problem of his famous 1900 
Paris address, D. Hilbert explicitly mentions Boltzmann’s limiting processes. But math¬ 
ematical progress was slow. Around 1950 H. Grad [6] first proposed a definite limit in 
which the Boltzmann equation would become exact, hence Boltzmann-Grad limit. One 
considers a fixed box Act 3 containing N hard spheres of radius e and has to study the 
limit where N —» oo, £ —> 0 such that the mean free path |A|/£ 2 1V remains constant. In 
this limit the volume filled by hard spheres tends to 0 proportional to e, hence physically 
the fluid is at low density. C. Cercignani [4] noted that Grad’s limit fits very nicely with 
the structure of the BBGKY hierarchy for hard sphere correlation functionqj. With this 
insight, at the 1974 Battclle Rencontres, Oscar Lanford advanced a complete proof [9]: 
he iterates the BBGKY hierarchy, regarding extra collisions with a given cluster as per¬ 
turbation. The resulting series is estimated uniformly for times t with \t\ < T, where T 
is of the order of (l/5)-th of the mean free time. For such a short time span one can 
thus rely on the term-by-term convergence of the series. In the Boltzmann-Grad limit 
the series corresponds to the iteration of the Boltzmann equation regarding the collision 
operator as perturbation. In summary the Boltzmann equation is obtained as a limit 
when starting from Newton’s equations of motion. The restriction to short times is an 
unavoidable feature of the perturbative expansion. To extend the validity of the result to 
longer times remains as an, apparently very difficult, open problem. 

But, even accepting the restriction to a finite kinetic time, Lanford’s theorem leaves 
many issues without answer. Can smooth potentials be included? There are no explicit 
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error bounds. Viewed probabilistically the Boltzmann equation is a law of large numbers. 
In fact, consider the number density on the one-particle phase space for a single many- 
particle trajectory. This approximates the number density computed from the Boltzmann 
equation, with a probability very close to 1 relative to the initial non-equilibrium measure. 
Is there then a central limit type result relative to this law of large numbers? 

The book under review presents a proof of Lanford’s theorem with precise bounds 
allowing both an explicit rate of convergence for hard spheres, and an extension to the 
case of smooth potentials. In the latter direction a first step was achieved already in the 
1975 Berkeley Ph.D. thesis by F. King [8]. As first observed by H. Grad [7], one has 
to write down a version of the BBGKY hierarchy which uses completed collisions. The 
differential form used conventionally in physics text books is of no use when inserted into 
the perturbation series. The properly modified BBGKY hierarchy contains couplings to 
all higher order correlations. But King proved the same uniform bound as Lanford and 
also that the terms beyond next order vanish in the limit e —> 0. It was then tacitly 
assumed that the term-by-term convergence could be accomplished by inspection. As the 
book demonstrates, this is far from the truth. Also, the explicit estimates require a much 
more detailed analysis even in the hard sphere case. 

Besides providing useful background material, the bulk part of the book is concerned 
with solving an n-particle scattering problem for arbitrary but fixed n. To understand 
the nature of this problem, we briefly have to outline the structure of the perturbation 
series. 

The starting point is a Hamiltonian system consisting of N particles which interact 
through a pair potential, 
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where (xi(t) , Vi(t)) e x R d , i — 1,... ,N, denote the position and velocity of particle 
i at time t , $ is a given two-body potential with range 1, and the physical dimension is 
d = 2,3. Particles are indistinguishable, they move on straight lines with constant velocity 
and are deflected at a mutual distance e. The singular case of hard spheres (formally, 
<h(a;) = oo if |x| < 1 and $(x) = 0 if |x| > 1) corresponds to the instantaneous elastic 
reflection. 

We shall fix Ne d = 1. Writing Z'n = (zi,..., Zn), Zi = the N— particle 

probability distribution is a function f N = f N (t,Z N ), invariant under permutations of 
the particle labels, with evolution governed by the Liouvillc equation. In other words, 
/jv(t) is obtained by transporting the prescribed initial distribution /(? along the solutions 
of the above system of ODEs. Integrating over all the position and velocity variables 
but the Erst s, one defines the s-particle marginal fffl. What we are interested in is the 
probability to find an arbitrary particle at position x and with velocity v at time f, i.e. 
f^\t,x,v). For N large, its evolution is expected to be well described by the Boltzmann 
equation 

d t f + v ■ W x f = [ [ [f'f[ - ffi\ b(v - v u u) dvidui, 
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where is the unit sphere in R d , f = f(v), f = f(v'), f[ = f(v [) and /i = f(v i) with 
v' — v + uj ■ (vi — v)uj, v[ = Vi — u j ■ (vi — v)uj. The collision kernel b(w,co)/\w\ is the 
differential cross-section of the two-body scattering. 

With the understanding that conditions on the potential <f> still need to be specified, 
a mathematical statement can be given as follows. 

Theorem. Let f 0 : R 2d —> M + be a continuous probability density such that 
ll/o exp(/3|u| 2 /2)|| L°° < +oo for some (3 > 0. Let E £ (Z S ) be the s-particle energy, 

E,(Z,)= £ J] 4 

l<i<s l<i<k<s 

Consider an initial “asymptotically independent” distribution of particles such that, for 
some fi > 0, the s-particle marginals satisfy |/^’ 0 (Z S )| < e us e -PE £ (z s ) un if orm iy { n 
Z s ,s,N, and fff'° —> f® s as N —>■ oo, locally uniformly outside the diagonals of physical 
space, {xi = Xk}- Then, there is a time T > 0 such that for t G [0, T] 

lim f£\t) = f(t), 

N^-oo 
Ne d ~ 1 =1 

where f(t ) is the solution of the Boltzmann equation with initial datum fo- 

A crucial byproduct of the proof is that, in the same limit, f^\t) —> f(t)® s for any 
fixed s. In particular, the asymptotic factorization of the state propagates for short times. 
The book provides a proof of the result, covering simultaneously the hard sphere potential 
and a certain class of smooth, short range repulsive potentials. The convergence is weak in 
the sense that one integrates against test functions with respect to the velocity variables, 
for any Z s outside diagonals of the physical space. In the case of hard spheres, a rate 
0(e a ) is shown to propagate in time, for any a < (d — l)/(d + 1) and Lipschitz / 0 . 

The classical strategy resorts to the full set of evolution equations for the family of 
marginals (f^^iKsKN- the BBGKY hierarchy. The rigorous derivation of these equations 
and the short time control, uniform in N , of their solutions, are discussed by dealing 
separately with hard spheres and smooth potentials (respectively parts II and III of the 
book). One faces indeed rather orthogonal difficulties. In the hard sphere case, a nontrivial 
discussion is needed to conciliate an integro-differential system of equations with the 
singular character of the dynamics. The iV-particlc flow is defined only away from a set of 
Lebesgue measure zero, while the collision operators are defined by integrals on manifolds 
of codimension 1. In the proper weak formulation the hierarchy reads 

f) fC) I r f( s ) _ n A s +1) 

u tJ N ' t-'sj N — ^ s ^s+lJ ]\r , 

where C s = • V Xj with boundary conditions identifying pre- and post-collisional 

velocities, and 

C s , s +ifN +i) = (■ N-s)e d ~ l Y, [ iy-(v s+1 -Vi)G <itS+1> flf +1) (-,Xi+£v,v s+1 ) dv dv s+1 . 

i=1 Jsf-'xRd' 
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Here G<j iS+ i> is an eventually negligible constraint ensuring that particle s + 1, with po¬ 
sition Xi + Ev, keeps distance larger than e from all the other particles, except i. Repeated 
iteration of the time-integrated hierarchy produces Lanford’s perturbative expansion. We 
refer to for the derivation of the perturbative setting. 

In case of a smooth interaction potential the scattering is not instantaneous, but space- 
time delocalized and the equations change drastically. The transport operator C s contains 
the additional term — X!i<i 7 y< s (l/ e )^ < ^(( a h — x j)/ £ ) • V Vi . The collision operator involves 
also velocity derivatives and the particle s +1 runs over the interior of a sphere of radius e 
around ay. The iterated expansion becomes cumbersome, since binary collisions are mixed 
up with many-body collisions. To overcome this problem, following King, one introduces 
truncated marginals, 


fV(t, Z.) 


l R 2d(N- a ) 


fN(t, Z s , Z s + 1, . . • , £/v) 11 ~^-\xi—Xj\>e d^s+1 • ■ ■ dzjy . 

ie{i,- ,s} 
je{s+i,...,N} 


These are asymptotically equivalent to the Moreover, in the hierarchy satished by 
{fff)i<s<N, the collision operator has the same structure of the hard sphere C S)S+1 , plus 
corrections containing higher order marginals and describing interactions of order three, 
four, etc., which a priori can be shown to be negligible in the Boltzmann-Grad limit 
(“cluster estimates”). The resulting perturbation series has the form 

°° rt fti ftk ~i 

In = / / •■■ / H s (i — ti) C S)S+ i H S+1 (H — t 2 ) C-s+i,s+2 

k=0 Jo Jo Jo 

... H a+ k{tk)f ( N +k),0 dt k ...dt 1 +0{e) , 

where H* is the s-particle flow associated to $ and /^ +fc ' ) ’° are the initial data. The 
0(e)-e rror reduces to zero only when returning to the hard sphere case. 

If f(t) is a solution to the Boltzmann equation, then the family of tensor prod¬ 
ucts f^(t, Z s ) := fit, Zi)f(t, z 2 ) ■ ■ ■ f(t , z s ) satishes formally an analogous representation 
(Boltzmann hierarchical expansion), with the operators H s and C SyS+ \ replaced respec¬ 
tively by the free flow operator and the hierarchical Boltzmann operator. This discussion 
is included in parts II-III of the book, together with the above derivations and a priori 
estimates, the notion of quasi-independence assumed for the time-zero state, and a precise 
formulation of main results in terms of hierarchies, see Theorems 8 and 11. 

Once the absolute convergence of the series expansion for fff is established for small 
times, one is left with showing the term-by-term convergence to the Boltzmann series. 
This is the subject of part IV. To understand the nature of the problem one has to look 
in detail at the structure of the generic term of the expansion. The term k = 0 is simple 
transport over the backward s-particle dynamics. The higher order terms are integrals over 
the trajectories of a special, fictitious, backward dynamics (in the book called “pseudo¬ 
trajectories”), in which k particles are adjoined to the first s ones sequentially at times 
t\ > t 2 > • • • > tk- Each new particle enters the evolution by appearing in a collision 
configuration, i. e. at distance £ from one of the already existing particles, its “progenitor”. 
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The precise way to add particles is specified by the integration variables appearing in 
the definition of the collision operator C s _ s+ \. Between t s+q and t s+q+ i, the backward 
dynamics H s+(? is applied. If the velocities of the new particle and its progenitor at time 
t s+q are post-collisional, then the pair will scatter under the action of H S+(J . The initial 
datum /^’ +fc ' ) ’ 0 is integrated over the time-zero values of these special characteristics. As a 
consequence, the problem can be rephrased as a result of convergence in the proper space 
of trajectories. Namely, the dynamics associated to the (modified) BBGKY perturbation 
series has to converge to the one of the Boltzmann series. The main difference between 
the two is that in the Boltzmann flow the only scattering events retained are those which 
happen at the addition of new particles at the times t\, £ 2 , ■ ■ •, bfc- In particular, each pair 
of particles collides at most once. Conversely, in the BBGKY flow other collisions, usually 
referred to as recollisions, may occur. Indeed, the heart of the proof consists in showing 
that the recollisions have a negligible measure in the space of pseudo-trajectories. 

I 11 the bulk of the book under review methods are developed which provide a quan¬ 
titative estimate of the recollision set. One can see this as a scattering problem: the 
recollision described by an operator C refers to the event that one of the two involved 
particles ends up hitting a third particle of the flow in the backward dynamics. This 
defines “untypical” values of the integration variables in C, which can be geometrically 
characterized and estimated. Such an estimate is rendered delicate by the fact that it 
requires a certain stability. The collision C itself introduces in fact additional deviations 
from the Boltzmann dynamics, due to a non-zero collision time and the associated spatial 
displacement. Furthermore, one must be able to iterate the process from C s+q to C s+q+ 1 . 
These difficulties are dealt with by a careful characterization of the set of untypical tra¬ 
jectories and a clever use of the dispersivity of the free flow. The technique applies in the 
same way to hard spheres and to the considered class of smooth potentials, where it leads 
to a rate of convergence depending implicitly on <f>. 

The results presented are a highly sophisticated completion of a program started some 
time ago. In kinetic theory, the conventional lore is that the Boltzmann equation is 
applicable for “any” potential which decays faster than Coulomb. For potentials of infinite 
range, no matter how rapid a decay, the current methods fail because the total cross- 
section diverges. But even for a finite range and smooth potential there will be initial 
conditions where the two-body collision time is infinite and the differential cross-section 
could be unbounded. The potentials covered in the book, see Sect. 15.2., have to satisfy 
a convexity assumption, which guarantees a bounded differential cross-section away from 
grazing collisions. Subsequent studies on the subject followed ra. which use geometrical 
arguments to control singular cross-sections and long collision times, thereby enlarging 
the class of admissible potentials. A more refined quantitative analysis for the hard 
sphere system provides information on cumulants through estimates of many-recollision 
events m- Recollision trajectories have also to be controlled over very long kinetic times 
for deriving Brownian motion of a tracer particle immersed in a rarified hard sphere gas 
close to equilibrium [2], see also pj for diffusion in the low density Lorentz gas. 

The presentation of the book is tuned so to mostly address an audience of mathe¬ 
maticians working on partial differential equations. Boltzmann’s 1872 effort, to conciliate 
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microscopic time-reversible dynamics with increase of entropy and trend to equilibrium, 

still remains as a source of challenging mathematical problems. 
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